The role of final-state interactions in Dalitz plot studies 
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Dalitz plot studies for multi-hadron decays of heavy mesons are expected to become very 
important tools for precision investigations of CP violation. A thorough understanding of 
the hadronic final-state interactions is a prerequisite to achieve a highly sensitive, model- 
independent study of CP-violating phases in such processes. We illustrate the theoretical 
tools available, as well as still to be developed, from low-to-medium-energy hadron physics 
for this purpose, and the goals of the informal Les Nabis network studying these and related 
problems. 



1 CP-violation in Dalitz plots 

A precise study of final-state interactions is increasingly becoming of paramount impor- 
tance for our understanding of the most diverse aspects of particle decays involving 
hadrons. Final-state interactions can be of significance for various reasons: if they are 
strong, they can significantly enhance decay probabilities; they can significantly shape the 
decay probabilities, most prominently through the occurrence of resonances; besides reso- 
nances, also new and non-trivial analytic structures can occur, such as threshold or cusp ef- 
fects (for the prominent role cusp effects have played recently in studying pion-pion inter- 
actions, see Ref. [1] and references therein); and finally, they introduce strong or hadronic 
phases or imaginary parts, the existence of which is a prerequisite for the extraction of CP- 
violating phases in weak decays (see e.g. Ref. [2]). Dalitz plot studies of weak three-body 
decays of mesons with open heavy flavor (both D and B) are expected to acquire a key 
role in future precision investigations of CP violation, due to their much richer kinematic 
freedom compared to the (effective) two-body final states predominantly used to study 
CP violation at the B factories. In many cases, the branching fractions are significantly 
larger; furthermore, the resonance-rich environment of multi-meson final states may help 
to enlarge small CP phases in parts of the Dalitz plot, and differential observables may 
allow to obtain information on the operator structure that drives CP violation beyond the 
Standard Model, once it is observed. Since the results from the B factories have shown that 
the Cabibbo-Kobayashi-Maskawa theory [3] represents at least the dominant source of CP 
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violation, our long-term goal will be to find other sources of CP violation that contribute 
additional, smaller effects. For this purpose, clearly extremely accurate measurements and 
means of theoretical interpretation are required. Strong evidence for CP violation in three- 
body final states has already been reported for B ± — > K ± 7r T 7r ± [4], e.g. with a 3.7c signal 
in the effective Kp channel, while only negative results exist for D decays so far [5]. 

There are different possibilities how to analyze CP violation in Dalitz plots. One sugges- 
tion is a strictly model-independent extraction from the data directly [6, 7], e.g. using the 
significance [7] variable defined as 



where N(i) and N(i) denote the event numbers of CP-conjugate decay modes in a specific 
Dalitz plot bin i. CP violation can then be identified in a deviation from a purely Gaus- 
sian distribution in the significance plots. The significance method allows to study local 
asymmetries and requires no theoretical input at all. 

An alternative approach, in contrast, makes use of information on the strong amplitudes 
as input. To see why this may be advantageous, we consider the following toy model: 2 
consider event number distributions given by a (Breit-Wigner) resonance signal (of mass 
M res and width r res ) on a certain background, with a CP-violating phase 5qp, according to 

(2) W ,N^ + ,Rei - } N — N — ^2,M„T„ 



M r 2 es + zM res r res J (s - M r 2 es ) 2 + (M res r res ) 2 • 

Figure 1 shows the count-rate difference and the corresponding significance distributions 
for two simulated cases of pseudo-data of different statistical weight: while in the high- 
statistics case, there indeed seems to be a deviation from Gaussian distribution in the sig- 
nificance, this is definitely lost in the case of less events. In contrast, fitting the data with 
the functional form (2) still allows to extract ^cp with some (limited) accuracy even in the 
sparser sample. So while the theoretical assumptions and prejudices going into such an 
analysis clearly have to be very carefully judged, their benefit in terms of vastly increased 
sensitivity is also obvious. 

In the following, we will therefore briefly sketch some of the tools available to analyze the 
hadronic amplitudes of the (light) final-state particles (such as pions and kaons). 



2 Pion-pion scattering: Roy equations 

Analyticity, unitarity and crossing symmetry provide a high degree of constraint for the 
pion-pion scattering amplitude, which can be exploited using dispersion relations. Start- 
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Figure 1 : Toy example for extraction of a CP-violating phase, for high- (left) and low- 
(right) statistics samples, using the known Breit-Wigner shape (top) versus the signif- 
icance variable (bottom). Input parameters used are Sqp = 5°, M res = 0.77 GeV, 
r res = 0.15 GeV. Figure courtesy of C. Hanhart. 



ing from a twice-subtracted dispersion relation at fixed Mandelstam variable t, 



(3) 



oo 

T(s,t)=c(t) + ± J ds'{* 



4Ml 



s' z {s'-s) s' z {s'-u) 



}lmT(s',t), 



the subtraction function c(f) can be determined from crossing symmetry. Projecting onto 
partial waves tj of definite angular momentum / and isospin I, one obtains a coupled 
system of partial-wave integral equations, 



2 oo . 

t I j(s)=k I J (s)+£ £ / ds'Kfj'^s^lmt^s'), 



(4) 



where the kernels Kj^s, s') are kinematical functions known analytically. The subtrac- 
tion polynomial fcj(s) contains the rcn scattering lengths as the only free parameters; these 
may in turn be further constrained by matching to chiral perturbation theory [8]. Equa- 
tion (4) can finally be turned into a coupled set of equations for the phase shifts, the Roy 
equations [9], by assuming elastic unitarity in the form 



(5) 



t}(s) 



, LmtJ(s) = a\t){s)\ 2 , cr=Wl 



4M2 
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Figure 2: Graphical representation of the consequence of analyticity and unitarity for 
form factors. 



Modern precision analyses of the Roy equations have been performed in Ref. [10], and 
similarly for pion-kaon scattering [11]; compare also the discussion of 77 — > tctc at this 
conference [12]. These provide us with high-precision parameterizations of the most rele- 
vant scattering amplitudes for light mesons appearing in the final states of heavy-meson 
decays. 



3 Analyticity and unitarity for form factors 

Final-state interactions between only two strongly interacting particles can be described in 
terms of form factors, which in turn can be linked to the properties of scattering amplitudes 
using analyticity and unitarity. As illustrated in Fig. 2, the unitarity relation for a form 
factor FJ(s) (here: of the pion) reads 

(6) Imf/(s) = F/(s) x 6(s - 4A4 ) x sin£j(s)<rW , 

from which one immediately deduces Watson's final-state theorem [13]: the form factor 
shares the phase 5j(s) of the (elastic) scattering amplitude. The solution to Eq. (6) is ob- 
tained in terms of the Omnes function [14], 

(7) F/(s) = P/(s)0}( S ) , Cr}(s) = exp J ^p^y} , 

where Pj(s) is a polynomial. Note that the Omnes function is completely given in terms of 
the phase shift. A classic application of such a form factor representation is the pion vector 
form factor Fy(s), which rather than by Eq. (7) is written in a more refined way as 

(8) F£(s) =n}(s)G w (s)n inel (s), 

where G a ,(s) takes into account p — <x> mixing, and Oinei( s ) parameterizes inelasticities, 
effectively setting in above s > (M n + M a ,) 2 , using conformal mapping techniques [15]. 
Such form factor representations can be used for analyses of the e + e~ — > n + n~ data to re- 
duce the error of the hadronic contribution to the muon g — 2, or to check the compatibility 
of the data with analyticity and unitarity [16]. Note that the effects of chiral dynamics are 
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Figure 3: Example for the complication of the analytic structure of 4-point functions 
through crossed-channel effects, here for the decay of an rj (double line) into three pions 
(single lines). 



particularly important for scalar form factors, where a parameterization in terms of Breit- 
Wigner resonances can lead to completely wrong phase motions (see e.g. Ref. [17] for the 
context of B — > 3n decays). 



4 Dispersion relations for three-body decays 

The application of dispersion relations to three-body decays is more complicated than the 
treatment of form factors due to the more involved analytic structure, and the possibility 
of crossed-channel rescattering (compare also Ref. [18] reported at this conference); see 
Fig. 3 for a depiction of the complication of the unitarity relation. We discuss here the 
(low-energy) example of r/ — > 3n decays, which has received much renewed attention 
recently [19,20] due to its importance for the extraction of the light quark mass ratios. One 
starts by decomposing the amplitude A4 (s, t,u) oc A(t] — > rc + n~ n°) into partial waves of 
isospin I according to [21,22] 

(9) M(s, t,u) = M {s) + (s - t)Mx(u) + (s - u)Mi(t) + M 2 (t) + M 2 (u) - ^M 2 (s) , 

where the A4i(s) are functions of one variable only, with only a right-hand cut. Equa- 
tion (9) is exact as long as discontinuities of D- and higher partial waves are neglected. 
The unitarity relation for the Mi(s), 

(10) ImA1;(s) = {Mifsj+^s)} x6>(s-4A4) x sin^(s)e-^ s ) 

(we now ignore the angular-momentum indices), is then complicated compared to Eq. (6) 
by inhomogeneities Aii{s), which are given by angular averages over the A4 j according to 

2 20 2 

M (s) = -(M )(s) + — (M 2 )(s) +2(s - s )(Xi)(s) + -k(s){zM 1 )(s) , 

l 

(z'7)(s) = \j dzz n f(l(3s -s + zk(s))) , s = i(Mj + 3M 2 n ) , 

-l 

(11) k(s) = y/(s- (Mr, + M n ) 2 ) (s - (M, - M n ) 2 ) x Jl-iM, 
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Figure 4: Real and imaginary part of the t] — > tt + tz tz° amplitude A4(s,t = u). Compare 
also Ref. [19]. 



and similarly for the other Kii. Note that the angular integration including the k(s) func- 
tion is non-trivial and generates a complex analytic structure, including three-particle cuts 
due to the fact that the rj is unstable and decays. The analog to the Omnes solution (7) are 
then integral equations involving the inhomogeneities [22] 

00 ^ 

/.^ k j / \ ^ /\f n 2 s 3 f ds' smS (s')Mo(s') ] 

(12) M,(.)-Q.(.)(«. + ». + 7.^ + - J ^ m ^J_t,e) }- 

with subtraction constants olq, fio, 70, and similar forms for the other partial waves. (See 
Ref. [23] for earlier, related formulations.) Equations (11) and (12) can then be solved itera- 
tively, e.g. matching the subtraction constants to chiral perturbation theory, see Fig. 4. The 
iteration converges fast, with the second iteration already very close to the final result. 

The method sketched here is currently also applied to other light-meson decays such as 
rf — > TjTCTt or co, cp — > 3n [24]. Challenges to be faced when extending this formalism 
to heavy-meson decays include the necessity to treat systems of integral equations when 
coupled channels within one partial wave cannot be ignored, or inelasticities are not neg- 
ligible. In particular when considering B-meson decays, elastic unitarity will surely not 
be sufficient. It will have to be checked when a purely perturbative treatment of crossed- 
channel effects is feasible (compare e.g. Ref. [25]), and when higher partial waves become 
important. To investigate these and related questions is part of the program of the infor- 
mal Les Nabis network [26], which brings together physicists from theory and experiment 
in heavy- and light-quark physics and aims at optimizing future Dalitz plot studies along 
the lines sketched here — with the strong goal to better interpret the mechanism of CP vi- 
olation in nature, yet at the same time teaching us important lessons on nonperturbative 
strong interactions. 
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